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Abstract 

In a bi-directional relay channel, two nodes wish to exchange independent messages over a shared 
wireless half-duplex channel with the help of a relay. In this paper, we derive achievable rate regions 
for four new half-duplex protocols and compare these to four existing half-duplex protocols and outer 
bounds. In time, our protocols consist of either two or three phases. In the two phase protocols, both 
users simultaneously transmit during the first phase and the relay alone transmits during the second 
phase, while in the three phase protocol the two users sequentially transmit followed by a transmission 
from the relay. The relay may forward information in one of four manners; we outline existing Amplify 
and Forward (AF), Decode and Forward (DF) and Compress and Forward (CF) relaying schemes and 
introduce the novel Mixed Forward scheme. The latter is a combination of CF in one direction and DF 
in the other. We derive achievable rate regions for the CF and Mixed relaying schemes for the two and 
three phase protocols. In the last part of this work we provide a comprehensive treatment of 8 possible 
half-duplex bi-directional relaying protocols in Gaussian noise, obtaining their respective achievable rate 
regions, outer bounds, and their relative performance under different SNR and relay geometries. 
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I. Introduction 

Bi-directional relay channels, or wireless channels in which two nodes (a and bj^ wish to exchange 
independent messages with the help of a third relay node r, are both of fundamental and practical interest. 
Such channels may be relevant to ad hoc networks as well as to networks with a centralized controller 
through which all messages must pass. From an information theoretic perspective, an understanding 
of these fundamental bi-directional channels would bring us closer to a coherent picture of multi-user 
information theory. To this end, we study bi-directional relay channels with the goal of determining 
spectrally efficient achievable rate regions and tight outer bounds to the capacity region. 

This two-way channel [4] was first considered in [22] where an achievable rate region and an outer 
bound for the case in which nodes operate in full-duplex were obtained. Non-orthogonaj^ full-duplex 
operation requires nodes to transmit and receive on the same antenna and frequency simultaneously. 
However, it may not be practically feasible to do so since the intensity of the near field of the transmitted 
signal is much higher than that of the far field of the received signal. In this work, we thus consider 
half-duplex communication in which a node may either transmit or receive at some time, but not both. 
Our goal is to determine spectrally efficient (measured in bits per channel use) transmission schemes 
and outer bounds for the half-duplex bi-directional relay channel and to compare their performance in 
a number of scenarios. These scenarios highlight the fact that different protocols may be optimal under 
different channel conditions. 

An obvious half-duplex bi-directional relay protocol is the four phase protocol, a — > r, r — > b, b— > r 
and r — > a, where the phases are listed chronologically. However, this protocol is spectrally inefficient and 
does not take full advantage of the broadcast nature of the wireless channel. One way to take advantage 
of the shared wireless medium would be to combine the second and the fourth phases into a single 
broadcast transmission by using, for example, network coding [1]. That is, if the relay r can decode the 
messages w a and Wb from nodes a and b respectively, it is sufficient for the relay r to broadcast w a © 
to both a and b. 

In this paper we consider two possible bi-directional relay protocols which differ in their number 
of phases. Throughout this work, phases will denote temporal phases, or durations. The three phase 
protocol is called the Time Division Broadcast (TDBC) protocol, while the two phase protocol is called 

'We call the nodes a and b terminal and source nodes interchangeably. 

2 By non-orthogonal we mean that no additional space, time, frequency, or coding dimensions are used to separate sent and 
receive signals. 



May 15, 2009 



DRAFT 



3 



the Multiple Access Broadcast (MABC) protocol. One of the main conceptual differences between these 
two protocols is the possibility of side-information in the TDBC protocol but not in the MABC protocol. 
By side-information we mean information obtained from the wireless channel in a particular phase which 
may be combined with information obtained in different stages to potentially improve decoding or increase 
transmission rates. The two considered protocols may be described as: 

1) TDBC protocol: this consists of the three phases a — ► r, b — > r and a <— r — > b. In this protocol, 
only a single node is transmitting at any given point in time. Therefore, by the broadcast nature of 
the wireless channel, the non-transmitting nodes may listen in and obtain "side information" about 
the transmissions of the other nodes. This may be used for more efficient decoding, i.e. improved 
rates. 

2) MABC protocol: this protocol combines the first two phases of the TDBC protocol and consists 
of the two phases a — > r <— b and a <— r — > b. Due to the half-duplex assumption, during phase 
1 both source nodes are transmitting and thus cannot obtain any "side information" regarding the 
other nodes' transmission. It may nonetheless be spectrally efficient since it has less phases than 
the TDBC protocol and may take advantage of the multiple-access channel in phase 1. 

We consider restricted protocols in the sense that the receivers must decode their messages at the end 
of the third phase (TDBC) or second phase (MABC) and collaboration accross multiple successive 
runs of the protocols are not possible. For each of the MABC and TDBC protocols, the relay may 
process and forward the received signals differently. These different forwarding schemes are motivated by 
different relaying capabilities or assumptions (about the required complexity or knowledge). Combining 
the relaying schemes with the temporal protocols, we can obtain various protocols whose rate regions 
are not in general subsets of one another. The relative benefits and merits of the two protocols and four 
relaying schemes are summarized in Tables U and JI] The four relaying schemes we consider are: 

1) Amplify and Forward (AF): the relay r constructs its symbol by symbol replication of the received 
symbol. The AF scheme does not require any computation for relaying, and carries the noise 
incurred in the first stage(s) forward during the latter relaying stage. 

2) Decode and Forward (DF): the relay decodes both messages from nodes a and b before re-encoding 
them for transmission. The DF scheme requires the full codebooks of both a and b and a large 
amount of computation at the relay r. 

3) Compress and Forward (CF): the relay does not decode the messages of a and b, nor does it simply 
amplify the received signal, but it performs something in between these two extremes. It compresses 
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the received signal, which it then transmits. To do so, the relay does not require the codebooks of 
the source nodes, but it does require the channel output distribution p(y r ) at the relay. 
4) Mixed Forward: the relay decodes and forwards (DF) the data traveling in one direction (from 
a — » b), while it compresses and forwards (CF) the data traveling in the opposite direction (from 
a <— b). For the mixed scheme, one of the codebooks and the channel output distribution are needed 
at the relay. 

In the CF scheme the relay searches the compression codebook to find an appropriate codeword. While 
the search operation is similar to the decoding operation in the DF scheme, the CF scheme may be less 
complex since the relay can choose a codebook for compression whose search space is smaller than the 
DF codebook. 

Some of these protocols and relaying schemes have been considered in the past. In [12], the DF 
TDBC protocol is considered. There, network coding in Z| is used to encode the message of relay r 
from the estimated messages w 3 and w\>. The works of [18] and [19] consider the MABC protocol, 
where an amplification and denoising relaying scheme are introduced. In [14] a lattice code is used for 
the Gaussian channel in the MABC protocol. The capacity region of the broadcast phase in the MABC 
protocol assuming the relay has both messages w 3 and u>b is found in [16]. In [23], [15] Slepian-Wolf 
coding is extended to lossy broadcast channels with side information at the receivers. In [10], achievable 
rate regions and outer bounds of the MABC protocol and the TDBC protocol with the DF relaying 
scheme are derived. There, network coding and random binning are the techniques employed to determine 
achievable rate regions. Uni-directional CF relaying in the full-duplex channel is first introduced in [5]. 
An achievable region in the CF MABC protocol is derived in [21]. In [8] a comparison between DF and 
CF schemes in full-duplex channels is performed, while in [20] a comparison of AF and DF schemes 
with two relays in the MABC protocol is performed. 

In this paper, we derive achievable regions for new CF and mixed relaying schemes in both the TDBC 
and MABC half-duplex protocols. We also obtain outer bounds for the TDBC and MABC protocols based 
on cut-set bounds. We compare the achievable rate regions of these four novel schemes with the regions 
and outer bounds derived in [10] as well as a simple AF scheme in Gaussian noise. We thus present a 
comprehensive overview of the bi-directional relay channel which highlights the relative performance and 
tradeoffs of the different schemes under different channel conditions and relay processing capabilities. 
Notably, we find that under some channel conditions the mixed TDBC protocol outperforms the other 
protocols and similarly, there are channel conditions for which the CF TDBC protocol has the best 
performance. 



May 15, 2009 



DRAFT 



5 



TABLE I 

Comparison between two protocols 



Protocol 


Side information 


Number of phases 


Interference 


MABC 


not present 


2 


present 


TDBC 


present 


3 


not present 



TABLE II 

Comparison between four relaying schemes 



Relaying 


Complexity 


Noise at relay 


Relay needs 


AF 


very low 


carried plus noise at rx 


nothing 


DF 


high 


perfectly eliminated 


full codebooks 


CF 


low 


carried plus distortion 


P(Vr) 


Mixed 


moderate 


partially carried 


one codebook, p(y,) 



This paper is structured as follows: in Section [TTJ we introduce our notation, review previously deter- 
mined achievable rate regions and outer bounds and define the protocols that we will consider. In Section 
ILTTI we derive achievable rate regions for the CF and mixed relaying schemes. In Section JV] we obtain 
explicit expressions for these, and previous rate regions and outer bounds in Gaussian noise. In Section 
Ivl we numerically compute these bounds in the Gaussian noise channel and compare the results for 
different powers and channel conditions. 

II. Preliminaries 

In this work, we will determine and compare the rate regions of eight bi-directional relay protocols. We 
consider the 2 phase Multiple Access and Broadcast (MABC) and the 3 phase Time Division Broadcast 
(TDBC) protocol versions of Amplify and Forward (AF), Decode and Forward (DF), Compress and 
Forward (CF) as well as a Mixed scheme which combines Decode and Forward in one direction with 
Compress and Forward in the other. The AF, DF protocol regions and the CF MABC protocol region 
have been derived in prior work [10], [19], [21] while the CF and Mixed protocol regions described in 
Section ITTT1 are determined here. Also we slightly improve upon the CF MABC protocol region in [21]. 
We formally define our notation and problem next. 



May 15, 2009 



DRAFT 



6 



TABLE III 



Input and output alphabets 



MABC 



TDBC 



Phase 1 



<-ta, x b , y r , y r 



x a ,y b ,y r ,y r l 



Phase 2 



Xr,y=,y b 



x b ,y a ,y„y. 



Phase 3 



N/A 



x r ,y 3 ,y b 



A. Notation and Definitions 

We consider two terminal nodes a and b, and one relay node r. Terminal node a (resp. b) has its own 
message that it wishes to send to the opposite terminal node, node b (resp. a). The relay node r may 
assist in the bi-directional endeavor. This paper will determine achievable rate regions for bi-directional 
relay protocols over half-duplex, discrete-time memoryless channels. The half-duplex constraint implies 
that a node cannot simultaneously transmit and receive data. 

We first start with a somewhat more general formulation of the problem that simplifies the application 
of cut-set outer bounds and then apply it the MABC and TDBC protocols considered here. We consider 
an m node set, denoted as M. := {1, 2, • • • , m} (where := means defined as). We use Rij to denote the 
transmitted data rate of message Wij from node i G M. to node j € M., i.e., Wij € {0, ... , [2 nRi - J J — 
1} := Si j. The protocols considered have either L = 2 (MABC) or L = 3 (TDBC) phases. We denote 
by > the relative time duration of the I th phase, where A^ = 1. For a given block size n, A^ n 
denotes the duration of the I th phase. Obviously, A^ n — > A# as n — ► oo. 

For notational convenience, we define the messages W 3 := W a ,b> := Wb,a an d the corresponding 
rates R 3 := i? a t> and := Rb,a- The two distinct messages W a and Wb axe, taken to be independent and 
uniformly distributed in the set of {0, ... , [2 nR '\ - 1} := S a and {0, . . . , [2 nRb \ - 1} := <S b , respectively. 

We use channel input alphabet X\ and channel output alphabet 3^ for node i. We will be constructing 
Compress and Forward schemes in which received signals are compressed or quantized before being 
re-transmitted. We let Y{ denote the compressed representation of the received signal at node i, which 
lies in the corresponding compression alphabet y% for node i. $?i is not necessarily equal to 3^. We 
summarize the input/output alphabets of the MABC and TDBC protocols in Table HlD In Section ITVl and 
|VJ we consider the case Xi = 3j = 3i = C, Vi 

For a given block length n, it will be convenient to denote the transmission at time 1 < k < n at node 

3 ^P r is used in the CF TDBC protocol only. 
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i by Xf, the reception at node i at time k by Y^. Note that the distributions of Xf and depend on 
the value of k, e.g. for k < Ai n • re we are in phase 1, for Ai n n < k < (Ai in + A2, n )^ we are in 
phase 2 and for (Ai n + A2 t n)n < k < n we are in phase 3 (in TDBC protocols only). During phase 
I we use x\ to denote the random variable with alphabet Xi and input distribution p^\xi). It is also 
convenient to denote by X$ := |i € S}, the set of transmissions by all nodes in the set S at time 
k, and by Xg := {X- ; |£ € S}, a set of random variables with channel input distribution p^\xs) for 
phase I, where xs '■= € S}. Lower case letters X{ will denote instances of the upper case Xi which 
lie in the calligraphic alphabets X^ Boldface Xj represents a vector indexed by time at node i. Finally, 
we denote xg := {xj|z € S} as a set of vectors indexed by time. In this work, Q will denote a discrete 
time-sharing random variable with distribution p(q). 

In order to define bi-directional communication rates we must define the encoders, decoders and 
associated probability of errors. We define Ws,t '■= {Wij\i € S, j € T, S, T C M}. For a block 
length n, let denote the set of received signals at node i up until the end of phase £, i.e., until 
the end of time n ■ Y^m=\ ^m,n- Encoders and decoders are functions xf^ ( W/ji m , y"^ ) G Xi and 
^Yj, , Wux^m) respectively for £ £ {1,2,- •• ,L} and for our proposed schemes, we will obtain 
single letter bounds. We define error events E^j := {Wij ^ ^ij(-)} f° r decoding the message Wij 
at node j at the end of the block of length re, and E\'- as the error event at node j in which node j 
attempts to decode wi at the end of phase I using a joint typicality decoder. 

Let A^\UV) represent the set of e-typical (uW,v") sequences of length n ■ A^ n according to the 
distributions U and V in phase The events flW(u,v) := {(u",v^) G A^(i7V)}. In general, joint 
typicality is non-transitive. However, by using strong joint-typicality, and the fact that for the distributions 
of interest x — > y — > y, we will be able to argue joint typicality between x and y by the Markov lemma 
of Lemma 4.1 in [2] and the extended Markov lemma (Lemma 3 of [17], Remark 30 of [11]). 

A set of rates Rij is said to be achievable for a protocol with phase durations {A^} if there exist 
encoders/decoders of block length n = 1, 2, . . . with both P[Eij] — > and A^ >n — > A^ as re — > oo for all 
£ An achievable rate region (resp. capacity region) is the closure of a set of (resp. all) achievable rate 
tuples for fixed {A^}. 

B. Previous results 

We use the following outer bounds and achievable rate regions of decode and forward protocols, derived 
in [10] for comparison purposes in Sections ITVl and IVl We simply state the results here for completeness. 
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Theorem 1: (Outer bound) The capacity region of the bi-directional relay channel constrained to the 
MABC protocol is outer bounded by the union of 

R a < mm{A 1 I(xi 1) ;Y r ^\X^,Q),A 2 I(xP;Y^ ) \Q)} (1) 
R b < mm{A 1 I(xl 1) ;Y r {1) \xi 1) ,Q),A 2 I(xP;Y a {2) \Q)} (2) 

over all joint distributions p(q)p^ (x a \q)p^ (x b \q)p^ (x r \q) with \Q\ < 4 over the alphabet X a x X b x X r . 

□ 

Theorem 2: (Outer bound) The capacity region of the half-duplex bi-directional relay channel con- 



strained to the TDBC protocol is outer bounded by 

R a < min^/CxiV^W 1 ^ (3) 

R b < m:m{A 2 I{X^M i \YP\Q),A 2 I{X { h 2) -M 2 \Q) + ^I{X^ ) -M :i) \Q)} (4) 

R 3 + R b < AiJ(X a (1) ; y r (1) |Q) + A 2 /(4 2) ; Y r (2) |Q) (5) 



over all joint distributions p(q)p ( - 1 \x a \ q)p ( - 2 \x b \ q) p^\x r \q) with \ Q\ < 5 over the alphabet X a x X b x X r . 

□ 

Theorem 3: An achievable rate region for the half-duplex bi-directional relay channel with the MABC 



protocol is the closure of the set of all points (R a ,R b ) satisfying 

R a < mm{A 1 I(X^;Y r {1) \X^,Q),A 2 I(X^;Y b {2) \Q)} (6) 

R b < min{Ai/(xJ 1) ; y r {1) |X a (1) ,Q),A 2 /(X r {2) ; y a {2) |Q)} (7) 

R a + R b < A 1 /(X a {1) , Xg) ; Y {1) \Q) (8) 



over all joint distributions p{q)p <yl \x a \q)p <yl \x b \q)p^ {x r \q) with \Q\ < 5overthe alphabet X a x X b x X r . 

□ 

Theorem 4: An achievable rate region for the half-duplex bi-directional relay channel with the TDBC 
protocol is the closure of the set of all points (R a ,R b ) satisfying 

R a <mm{AJ(xi 1) ;Y r {1) \Q),A 1 I(xi 1) ] Y b {1) \Q) + A 3 I(X^ (9) 
R b < min { A 2 I(X {2) ; Y r {2) \Q), A 2 I{X {2) ; Y a {2) \Q) + A 3 /(X r (3) ; Y a {3) \Q) } (10) 

over all joint distributions p(q)p^' (x a \q)p^ (x b \q)p^ (x r \q) with \Q\ < 4 over the alphabet X a xX b x X r . 

□ 
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C. Compress and Forward using two joint typicality decoders 

In Compress and Forward protocols, unlike in Decode and Forward protocols, the relay node r does 
not decode the message w 3 or Wb. Thus, network coding techniques such as the algebraic group operation 
w 3 Wb used in [10] cannot be used to generate w r for the current CF schemes. Instead, two jointly 
typical decoders at each node are used to decode w r . 




Fig. 1. The data flow in the compress and forward MABC protocol 

To illustrate the decoding scheme, consider the decoder at node a which wishes to decode the relay 
message w r in order to ultimately decode the desired message from node b, Wb- After phase 2, node a has 
the known sequences x^^Wa) and y a . Node a then finds the sets of all y^(u> r ) and xi 2 \w r ) such that 
(xi 1 ^ (w a ), fr (w r )) and (~x.i 2 \w r ),y^) are two pairs of jointly typical sequences, as shown in Fig. Q] 
Then node a decodes w r correctly if there exists a unique w r such that (x^ {w a ), Yr (w r )) S A^ 1 \X a Y r ) 
and (xp-^uvXy^) e A^(X r Y a ) and declares a decoding error otherwise. 

III. Achievable rate regions for Compress and Forward and Mixed Protocols 

In this Section we present three new achievable rate regions in Theorems [8j [12] and Q31 and a slight 
improvement of [21] in Theorem [5] Theorems [5] and [8] are for two phase Multiple Access and Broadcast 
(MABC) protocols. In MABC protocols, nodes a and b transmit simultaneously as in a standard multiple- 
access channel in phase 1, the relay processes the received signal (either by decoding, amplifying, or 
compressing the signal as dictated by the protocol), and during phase 2 the relay broadcasts its signal to 
the two nodes. In Theorem [5] the relay simply uses a CF operation for both messages/directions, while in 
Theorem [8] the relay uses DF to transmit the w a message while it uses CF to transmit the Wb message. The 
final two theorems [12] and [14] employ three phase Time Division Broadcast (TDBC) protocols. During 
the first phase, node a transmits while both the relay and node b receive its signal. During phase 2, 
node b is the sole transmitter while node a and the relay receive its transmission. After these phases, 
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a 


r 


b 


phase 1 






• 


• 









compress : y? 1 — • yrJWr) 
encode : Wr — • xf 2 ' 




phase 2 









• - 


• 





decode : xS! yf 1 — »Wr i 


! decode 


xfi! yS^-Wr 


decode : x§! yP 1 — <wb 1 


i decode 


xfi! yPUwa 



Fig. 2. The two-phase MABC protocol with a relay using a CF scheme. 



the relay processes the received signals and is the sole transmitter during phase 3, during which it can 
aid nodes a and b to recover each others' messages. We now proceed to describe each protocol more 
precisely and present their respective achievable rate regions. Again, in Theorem [12] the relay uses CF 
when re-transmitting both messages, while in Theorem [141 the relay uses DF to transmit message w a and 
CF to transmit message w\>. We now proceed to the main technical results of this work. 

A. MABC Protocol 

In the MABC Protocol, message w a is communicated from node a to node b and message Wb is 
communicated from node b to node a with the help of the relay in two phases as shown in Fig. [2] and 
[3] During phase 1, nodes a and b simultaneously send independent messages w a and Wb as codewords 
**\w a ) and x^ 1 (i^b) to the relay, forming a classical multiple-access channel. Since we assume half- 
duplex nodes, neither a nor b can receive the message of the other during phase 1. The relay receives 
the signal y^ according to p{yr\xa,x9'). Rather than attempting to decode message w a and (as 
in a DF scheme), it compresses the received yP into a signal yf 1 ' (w r ). The index w r is then mapped 

(2) 

in a one-to-one fashion to the codeword x|: (w r ) which is broadcast in phase 2 back to the relays. The 
challenge here is to determine the optimal compression strategy such that just enough information is 
carried back to the nodes to decode the opposite node's message. A key observation is that the nodes 
may use their own phase 1 transmitted messages as side-information in the decoding of phase 2 signals. 

Theorem 5: An achievable rate region of the half-duplex bi-directional relay channel with the compress 
and forward MABC protocol is the closure of the set of all points (R a , R^) satisfying 

R a < AJixPiYPlX^^) (11) 
i2 b <AiJ(^ 1) ; y r (1) |X a (1) J Q) (12) 
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a 


r 


b 


phase 1 








• 





V 




decode : yr — 'Wa 
compress : y^— > y'lwro) 
encode : Wa.wro — .uf 21 

WrO >ut 21 

Wa, WrO >Xr 21 




phase 2 









• - 


• 





decode 


vin ,,( 2 >_ „ ! 

Xa, ya 'WrO I 


1 decode : yd 2 '— «wa 


decode 


Xa! yr' 1 — >Wb 1 





Fig. 3. The two-phase MABC protocol with a relay using a mixed scheme. 



subject to 

Ax/^y^lX^Q) < A 2 /(X r (2) ; y b (2) ) (13) 
Ax/^Y^lxi^Q) < A 2 I(xP;Y a i2) ) (14) 

over all joint distributions, 

p(q, x a , x b , x r , y a ,y b , y r , Vr) = P {1) (q, x a ,x b ,y n y r \q)p {2) {x r , y a , y b ) (15) 

where 

p (1) (q, x a ,x b ,y r , y r ) = {q)p W (x a \q)p W (x b |<?)p (1) (y r \x a , x b )p {l) {y r \y r , q) (16) 
P^ 2) (x r ,y a ,y b ) = p( 2) {x r )p( 2) (y a ,y b \x r ) (17) 

with \Q\ < 4 over the alphabet X a x X b x X r x y a x y b x y r x j) r . □ 
Remark 6: The bound of Theorem [5] is essentially derived in [21]; equation (fT6l ) is a slight extension, 
as we use p^(y r \yn q) instead of (y r \y r ), i.e., in [21] the codewords y^ are generated according to 
P^ {Vr) = YlP^iu^P^Hyrlyr), while in (fToT ) the distribution space p^ l \y r \q) = Yl {Vr)p^ {vAVr, q) 
is larger. By conditioning on q, one can "fine-tune" the distribution of for each given q and the left 
side of ( fl"3l ) and (fT4l) can be reduced. This is because the distributions of X a ^ and Xl , and hence 
Y"/ 1 ^, depend on For example, let p^(q = 1) = a n and p^(q = 2) = 1 — a„, where < a n < 
1. For q = 1 we optimize p^(yr|l) and generate (a n Ai jn • n)-length sequence yr (w r i), w r \ € 
{0, l,--- L2 ni?rl J}, where R rl = a n A ltn (I(Y r (1) ; Y r (1) \q = 1) + e). Likewise, we generate y r (1),2 K 2 ) 
for q = 2. To compress yr to y£ , we construct y^ = (y£ ,yr 2 ) and y^ = (yr ^yr ) and 
choose w r = (w r i, '^r2) if both (yr , yr (wri)) and (y^' 2 , yr 2 (i«r2)) are jointly typical. Then the 
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rate R r = R r i + R r2 = Ai/(y r (1) ;Y r (1) |g). However, if one generated y r from the distribution p^(y r ) 
then R r = AiI(Y T ^ ;Y T ^) > A\I(Y r ,Y T ^\Q) with strict inequality except in degenerate cases. 

Remark 7: Strong typicality is required for the proof of Theorem [5] in order to apply the Markov 
lemma to (X^jX^) — > y} 1 ^ — ► Y^ for each given q. Since strong typicality is defined for discrete 
alphabets, Theorem [5] cannot be directly extended to continuous alphabets. However, the extended Markov 
lemma (see Remark 30 of [11] as well as Lemma 3 of [17]) shows that for Gaussian distributions, the 
Markov lemma still applies. 

The previous theorem assumed the relay used a CF scheme for both messages w a and w b . However, in 
the event of asymmetric channel gains between the two nodes at the relay, it may be beneficial to have the 
stronger channel use a DF scheme while the weaker channel uses a CF scheme since decoding may not 
be possible. We next consider a mixed MABC strategy in which phase 1 is still a multiple access channel 
and phase 2 is still a broadcast channel. However, the relay uses a novel strategy in determining its phase 

(2) 

2 codeword Xr . That is, the relay operates such that the a — ► r — ► b link uses decode and forward while 
a <— r <— b link uses compress and forward. Furthermore, the relay applies a Gelfand-Pinsker coding 
scheme to protect w a in the r — > b link. In this case, an achievable rate region is given by Theorem [8] 

Theorem 8: An achievable rate region of the half-duplex bi-directional relay channel with the mixed 
forward MABC protocol, where a — > b link uses decode and forward and b — > a link uses compress and 
forward, is the closure of the set of all points (R a ,R b ) satisfying 

R a <mm{A 1 I(X< 1) ;Y r (1) |Q),A 2 /^ (18) 
R h < A 1 I(xl 1) ;Y r {1) \xP,Q) (19) 

subject to 

AiJ(y r (1) ;y r (1) |JT a (1) ,Q) < min{A 2 /([/ r (2) ,C/( 2) ; y a (2) |Q), A 2 /(^ 2) ;C/ r (2) ,y a {2) |Q)} (20) 
over all joint distributions, 

p(q, x a ,x h ,x r ,u a ,u b , u r , y a , y h , y r ,y r ) = p(q)p {1) ( x a i yri 2/r 

\q)p^(u b ,u r ,x r ,y a ,y b \q) (21) 



where 

p W {x a ,x b ,y r ,y r \q) = p W (x a \q)p W {x b \q)p W (y r \x a , x b )p (1 \y r \y r , q) (22) 

P {2) {u b , u r , x r , y a ,y b \q) = p {2) {u b , u r \q)pW (x r \u b , u r , q)p {2) \y a ,y b \x r ) (23) 

with | Q\ < 7 over the alphabet X a x X b x X r x U b x U r x y a x y b x y r x j) r . □ 
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Remark 9: In the second phase, the relay broadcasts the received signals from the first phase. In contrast 
to the DF and CF schemes, one of the terminal nodes (in this case, node a) has perfect information of 
interference at node a, while there remains unknown interference at the other side (at node b). We use a 
Gel'fand-Pinsker coding scheme [7], [6] for the link r — ► b which yields the second term of (TT8T ). From 
the side information w a available at node a, node a is able to reduce the interference, yielding (T20b . 

Remark 10: If we apply the achievable bound to the Gaussian noise channel without fading with 
Costa's setup in [3] with \Q\ = 1 we have in phase two: 

= Vr {2) + atf > (24) 

= F r (2) + Ut 2) + ZP (25) 

y b (2) = + + Z« (26) 

where Z {2) ~ CM(0,N a ), Z {2) ~ CM(0,N b ), U {2) ~ CM(0,P Ur ), V r {2) ~ CM(0,P Vr ), and V r (2) , U (2) 

(2) 

are independent and Pjj r + Py r = P r . V r is an intermediate random variable generated according to 
£A/"(0, P r ) which contains the information to be transmitted from r to a. Then from (fT~8T > and (|2Qb the 
achievable rate of link r — > a, i? ra (resp. i? r b for r — ► b) is : 

*ra = mm |log 2 [1 + - j , log 2 ^ _ j | (27) 

o I ( P Vr (P Vr + P Ub +N b ) \ 

Rrh ~ l0g2 {PvP Ub (l-*) 2 + N b (P Vr + aiP Ub )) (28) 
Proof: Random code generation: For simplicity of exposition, we take \ Q\ = 1. 

1) Phase 1: Generate random (n • Ai n )-length sequences 

. xi 1} (u/ a ) i.i.d. with p«(x a ), io a G <S a = {0, 1, • • • , [2 nR >\ - 1} 

. xf^twb) i.i.d. with pWfcb), ttfb G 5i, = {0, 1, ■ ■ • , L2 niJb J - 1} 



. y r (1) Ko) i-i-d. With p«(y r ) = E^P (1) (l/r)P (1) (yr|l/r) , Wrf) € {0, 1, • - - , [2 nR "> ] - 1} := 5, 



r0 



2) Phase 2: Generate random (n • A2, n )-lengfh sequences 

. u[ 2) (w r ) i.i.d. with p( 2 ) (u r ), tu r € {0, 1, • • • , [2 nR '\ - 1} := 5 r 
• u b 2) (-u; r o) i.i.d. with p (2) (ii b ), w r0 G S r0 
and define bin B, L := {w r \w r G [(» - 1) • [2 n ^-^)j [2 n (^-^)j]} for i G 5 a . 

Encoding: During phase 1, the encoders of node a and b send the codewords x a (u; a ) and x b ^(w;b) 
respectively. At the end of phase 1, relay r decodes w a and maps y£ to a message index w r o if there 
exists a w r Q such that (yp ,y^(^ r o)) £ -A^O^r)- Such a u) r0 exists with high probability if 

i? r0 = Ai, n /(y r (1) ; y r (1) ) + e (29) 
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and n is sufficiently large. We choose 

i2 r = A 2 , n /(£/ r (2) ;Y b (2) )-4e. (30) 

To choose w r , the relay first selects the bin Bw a and then it searches for the minimum w r G Bw, such 
that (ui 2 \w r ),u^ (w r o)) G (U r Ub). This ensures uniqueness of w r if such a w r exists, i.e., w r is a 
function of {w a ,w r o). Such a w r exists with high probability if 

|^-,J>2"'( A -^ (2, ^ <2, ) +e ). (31) 

Since \B{\ = 2 n ( R, ~ R *\ \/i G S a , this condition is equivalent to 

R a < A 2 , n /([/ r (2) ; y b (2) ) - A 2 , n /([/ r (2) ; U {2) ) - be. (32) 

The relay then sends randomly generated i.i.d. according to p^ 2 \x r \u r ,u\^) with u?'(w r ) and 
u b v^r-o) during phase 2. 

Decoding: Node a estimates iZ? r o after phase 2 using jointly typical decoding. First, since a knows 
w a , it can reduce the cardinality of w r to |^2 n (- Rr_ - Ra )j. Furthermore, it forms two sets of w r $ based 
on typical sequences, {w r0 \(yi a l \w a ),yr\w r o)) G ylW(X a y r )} and {w r0 \(ui 2) (w r ), u fa 2) (w r0 ), y a 2) ) G 
A^ 2 ) (C/ r f7by a )i G -B«, a }- After decoding w r o (which is a success if there is a single common element in 
both of the previous sets), a decodes w\> using jointly typical decoding of the sequence (xa , Xu ,yr ). 
Node b decodes w r after phase 2 and from the bin index of w r it estimates w a . 

Error analysis: By the union bound, 

P[£ a , b ] < P[E$ U E®] (33) 

<P[E$] + P[E®\Eg;)] (34) 

P[E b>a ] < P[E$ U E®] (35) 

<P[E$]+P[E®\E$] (36) 

We define error events in each phase as follows: 

1) E$=E$> 1 UE$'\ 

ES' 1 : (x«K),y«)^A«(X a y r ). 

Ei]}' 2 : there exists w a + w a such that (x£\w a ),yP) G A^\X a Y r ). 

2) ^^Jfu^U^U^. 

-E^b' 1 : there does not exist a w r0 such that (y r ,y^(w r0 )) G -A^OW)- 

i*Vb' 2 : there does not exist a ii; r G -B Wa such that (u^(to r ),uJ~(u;ro)) G -4 (2) (^r^b)- 
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E r ,a 

p(2),3 

P (2),4 
^r,a 

P (2),5 



E r (2) ' 4 : there exists w r / w r such that (u[ 2) (w r ), y b 2) ) G A( 2 )([/ r y b ). 

3) E® = fig?' 1 U Eg' 2 U Eg' 3 U Eg' 4 U Eg' 5 U Eg' 6 . 

Eg' 1 : there does not exist a w rQ such that (y^.y^Ko)) G A^(Y r Y r ). 

there does not exist a u; r G B Wa such that (ur (u; r ),u b (w r o)) G yl( 2 )([7 r i7 b ). 

(u?\w r ),u®(w r0 ),y { a 2) )?A( 2 Xu r u b Y a ). 

(2) f2l (2) 

there exists (w r ,w r o) where w r / w r andtZ; r o ^ uvo such that (uj (w r ),u b ; (uvo),y a ) G 
A( 2 )({/ r {7 b Ya) and (x^ (w a ),y r (1) (w r0 )) G AW(X a y r ). Recall, w r is uniquely specified by 
(w a , w r o). Hence for a given w a , there are at most 2 nRr0 such (w) r , i^o) pairs. 
Eg' 6 : there exists w r0 ^ w r0 such that (u r (2) K), u b 2) (w r0 ), y { a 2) ) G ^( 2 )([/ r C/ b y a ), 
(xi 1) K),y r (1) Ko))GA( 1 )(X a y r ). 

4) Eg = Eg' 1 U Eg' 2 - 

Eg' 1 : (x^^xkVbU^Ko)) ^ A«(X a X b y r ). 

Eg' 2 : there exists w b / w b such that (x^Va^x^^y^Ko)) G (^ a X b y r ). 

Then: 

P[E$) SPIES' 1 ] + P[Ei]}' 2 ] (37) 
=P[5W(x a K) Jr )] + P[U 4M J)W(x a («i a ),y r )] (38) 

< e + 2 ™(^-Ai,„/(Xi 1) ; y r <1, )+3e) (39) 

E[E; 2) |E«] ^EfEg' 1 !^] + P[E r (2 j' 2 |E«] + P[E^' 3 \E^] + P[E r (2 b ),4 |E«] (40) 

<2e + E[E( 2 )(u r K),y b )] +E[U^E( 2 )(u r K),y b )] (41) 

< 3e + 2 ™(ft-A 2 ,„/(C7 r (2) ; y b <2) )+3e) (42) 

In (HB, EfE^' 1 !^ 1 ^] and P[E r ^' 2 |E^r ] are less than e due to ([29]) and ([32]), respectively. Furthermore, 

P[Eg] ^P^g' 1 ] + P[Eg' 2 ] + P[Eg' 3 ] + P[Eg' 4 ] + P[Eg' 5 ] + P[Eg' 6 ] (43) 
<2e + P[P>W (x a (w a ), y r Ko))] + Fp (2) (u r ( Wr ),u b ( Wr o),y a )] + 
P[U ^ E( 2 )(u r (w) r ),u b ( ? Z; r o),y a ),E( 1 )(x a K),y r (u; r o)))] + 
P[U «,=„, E^ (u r K), u b (u; r0 ), y a ), E« (x a (w a ), yrCwrf,)))] (44) 



<4 e + 2™(fto-A 2 ,„/(C/ r (2) ,C/ b (2) ;y a <2) )-A 1 ,„/(y r (1) ;X a (1) )+7e) + 
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2 n( J R r0 -A 2 .„/(C/ b <2) ;l/ r <2) ,y a (2) )-A 1 ,„7-(y r (1) ;Xi 1) )+7e) ( 45 ) 

P[E$\E$) <P[E { £ A \E^} + P[E®> 2 \E$] (46) 
=P[fl (1) (x a ( Wa ),x b K),y r ( Wr0 ))] + J P[U^ u , bJ D (1) (x a K),x b (u; b ),y r (u» r0 ))] (47) 

< e + 2™(^-Ai, n /(X b <1) ;K <1) |X< 1) )+4 £ ) (48) 

In ([44]), P^J' 1 ] and P[P r ( , 2) ' 2 ] are less than e due to (|29]> and d32]), respectively. In (|45]>, PfP/W (x a (w a ), y r (u> r o))] 
is less than e by the Markov lemma and the total cardinality of the case (w r ^ w r and w r Q ^ w r o) is 
bounded by 2 nR, ° since w r is uniquely specified if (w a ,w r o) is given. 

Since e > is arbitrary, the conditions of Theorem [8j d30l and the AEP property guarantee that the 
right hand sides of ([39]), (02]), (05]) and (08) corresponding to the first term of (HD, d30]), (f20]) and (fl9l ) 
respectively vanish as n — > oo. By the Caratheodory theorem in [9], it is sufficient to restrict \Q\ < 7. ■ 

The mixed MABC region in Theorem [8] is outer bounded by the DF MABC region in Theorem [3] In 
the mixed MABC protocol, the relay r has to be able to decode w a correctly after phase 1 without any 
information about u? b . If node a can decode w; b from a compressed version of y r and knowledge of w a , 
then by the information processing inequality, node r can decode from y r and w a . 

Theorem 11: The achievable rate region of the mixed MABC protocol of Theorem [8] is outer bounded 
by the achievable rate region of the DF MABC protocol of Theorem [3] □ 
Proof: If R a and P b lies in the region of the mixed MABC protocol of Theorem [8] for a given 
distribution (x a (x b )p (1) (Ur \x a , x b )p^ (y r \y r )p^ (x r )p^ (y a ,yb\x r ) then: 

R a < A 1 I(xi 1) ;Y r {1) ) < A!/(X a (1) ;y r (1) |X«) (49) 
P a < A 2 /(C/ r (2) ; y b (2) ) - A 2 /(C/ r {2) ; U®) < A 2 /(X r (2) ; Y b (2) ) (50) 
P b < AxI(xW; Y r {1) \X^) < A^X^ y (1) |X a (1) ) (51) 

p b < wxfWW) < WCrPWW) 

< A 2 I(uP\u^;Y a {2) ) < A 2 I(xP;Y a (2) ) (52) 
R a + P b < A 1 I(XP;YP) + A^xW^^lxP) < AiJ(X a (1) , if 1 ; Y" r (1) ) (53) 

(ED and are from the Markov process (P r (2) , f/^ 2) ) -> X r (2) -> (Y a (2) , Y b (2) ). From (09]) - (El, P a 
and P b are in the region of the DF MABC protocol in the Theorem [3] Therefore, every point (P a ,P b ) 
represented by a convex combination of distributions (or using Q) in the Mixed MABC protocol is also 
achievable with the DF MABC protocol. ■ 
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From Theorem [TT] the mixed MABC protocol does not achieve any rate pairs which cannot be achieved 
by the DF MABC protocol. However, its possible benefit lies in that it only requires the relay r to possess 
one of the codebooks of a and b. Therefore, in the event that relay r has one of the codebooks of the 
terminal nodes, by employing the Mixed MABC protocol one can achieve a rate region which outperforms 
that of the CF MABC protocol. In practice, if there are many terminal and relay nodes a relay may have 
some but not all codebooks of the terminal nodes. If the relay has full codebook knowledge then DF is 
possible, otherwise CF (or AF) may be more appropriate. 

B. TDBC Protocol 

The Time Division Broadcast (TDBC) protocol consists of three phases rather than the two seen in 
the MABC protocol. The MABC protocol takes advantage of the various gains provided by multiple- 
access schemes by having both nodes a and b transmit during phase 1 . However, the possible direct links 
between nodes a and b are not exploited. The TDBC protocol aims to exploit the direct link by having 
the nodes combine the signals received on the direct link and through the relay node, that is, exploit the 
side-information available at the decoders. 

The TDBC protocol consists of three phases, as illustrated in Fig. [4] and [5] During phase 1, node a 
is the sole node to transmit, while both the relay and node b receive this transmission. During phase 
2, node b transmits while the relay and node a receive. After phase 2, the relay processes the signals 
received during the first 2 phases and proceeds to broadcast to nodes a and b during the third phase. 

In the CF TDBC protocol, we use two different broadcasting schemes in the last phase. For this reason, 
we divide phase 3 into two sub-phases. In the first relay-broadcast phase, we use Marton's broadcast 
scheme of [13], in which two different messages are transmitted to the two receivers. In this scheme, 
neither receiver uses side information (w a at node a and w\> at node b) to decode the messages. In the 
second relay-broadcasting phase, we assume a compound channel, i.e., a common message is transmitted 
to the two receivers which have different side information. 

For convenience of analysis, we denote the first part of the relay-broadcast phase as phase 3 and the 
second as phase 4. 

Theorem 12: An achievable rate region of the half-duplex bi-directional relay channel with the com- 
press and forward TDBC protocol is the closure of the set of all points (R a , i?t>) satisfying 



R a < Ai/(x a (1) ; y r (1) ,y b (1) |Q) 



(54) 



R b < A 2 /(xf ; y r (2) ,y a (2) |Q) 



(55) 



May 15, 2009 



DRAFT 



18 



phase 1 •- 



compress : — • yfJWao) ! 



phase 2 • 



compress : yr'— »yr 2 Wbo) ! 
encode : w a o,Wbo — *Ua'! 
encode : w a o,Wbo — tub 31 ' 



phase 3 



decode : yi^Wbi ! encode : w a o,Wbo — »xf 4 ' ! decode : yb <3L -»Wai 



phase 4 



decode : xS,'yf,'yS!wbi- 
decode : ya| y { r 2 '— »wb 



decode : xf^ybVyb'lwai— »wao 
decode : ytij y'r 1 '— *w a 



Fig. 4. The three-phase TDBC protocol with a relay using a CF scheme. 



a 


r 


b 


phase 1 








• 


• 







! decode : yr'-^wa ! 
■ ■ 


phase 2 


• 






• 











compress : y?— > yf 2 wro) 
encode : Wao. Wro — >uf 3 ' 

WrO • ui?> 

WaO, WrO 'Xp 1 




phase 3 










• - 


• 






decode 


Wa, ya', ya '— 'WrO 




decode 


yo 3 '— 'WaO 


decode 


.») .-.(2) ,7/u 
ya, yr — »Wb 




decode 


yS! WaO — >Wa 



Fig. 5. The three-phase TDBC protocol with a relay using a mixed scheme. 



subject to 

a a AiJ(F r (1) ;F r (1) |Q) < A 3 /([/ a (3) ;F b (3) |Q) (56) 
a b A 2 /(F r (2) ;F r (2) |Q) < A 3 /(^ 3) ;F a (3) |Q) (57) 

aaA^F^F^IQ) + a b A 2 /(F r (2) ;F r (2) |Q) < A 3 /([/ a (3 ^ F b {3) |Q) + A 3 /(^ 3 ^ F a (3) |Q) - A 3 /(f/ a (3) ; ^ 

(58) 

(1 -a a )A 1 /(F r (1) ;F r (1) |Q) + A 2 /(F r (2) ;F r (2) |4 2) ,Q) < A 4 I(X^ ;Y^) + A 1 I(Y r W ;Y^\Q) (59) 
(1 -a b )A 2 /(F r (2) ;F r (2) |Q) + A^F^; F r (1) |X a (1) , C?) < A 4 /(X r (4) ; F a (4) ) + A 2 /(F r (2) ; F a (2) |Q) (60) 
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where < a a , a b < 1 over all joint distributions, 

p(q,x a ,x b ,x r ,y 3 ,y b ,y n y r ) 

= p{q)p^ (x a ,y h ,y r , yM)v {2) {xb, y a ,y r , yM)p {3) {u a ,u b , x r , y a ,yb\q)p {i Hx r , y a , y b ) (61) 

where 

p (1 \x a ,yb,y r ,y r \q) = p (i: \x a \q)p^ l] '(y b ,y r ka)p (i; \y r \y r ,q) (62) 
P {2 \x b ,y a ,y r ,y r \q) = p i2) (x b \q)p^\y a ,y r \x b )p^ 2) (y r \y r ,q) (63) 
p (3) ( u a, u b , x r , y a ,yb\q) = p^\u a ,u b \q)p^\x r \u a ,u b , q)p {? \y a ,yb\x r ) (64) 
P (4) (^r,y a ,yb) = P^(x l -)p^(y a ,y b \x r ) (65) 

with \Q\ < 13 over the alphabet X a x X b x X? x y^ x x yf x j> r 2 . □ 
Remark 13: If side information is very limited, then with A4 — > 0, a a , a b — > 1 the relay phase acts as 
a classical broadcast channel. At the opposite extreme, the side information cancels out all interference, 
i.e., /(Y r (1) ;Y r (1) |X a (1) ) = I(Y r {2) ;Y r {2) \xl 2) ) = 0. Then we set A 3 -> 0, a a ,a b -> 0. H 
Proof: Random code generation: For simplicity of exposition, we take [Q| = 1. 

1) Phase 1: Generate random (n • Ai „)-length sequences 

. xi 1 ^) i-i-d- with pWfca), w a G «S a = {0, 1, - ■ ■ , [2 nR >} - 1} 

n -1 I onR.n I 11 . t 

>a0 

and generate a partition of S a o randomly by independently assigning every index w a o G 5 a o to a 
set S 3 q i, with a uniform distribution over the indices i G {0, . . . , L2 ni?al J - 1} := S al . We denote 
by Sao(^ao) the index i of 5 a cM to which w a o belongs. 

2) Phase 2: Generate random (n • A2, n )Tengfh sequences 
. x i2) (w b ) i-i-d. with p^(x b ), w h G «S b = {0, 1, ■ • • , L2 n ^J - 1} 

• yr (2) Ko) i-i-d. with pM(y r ) = P {2) {y r )p {2) {$r\Vr) , w M G {0, 1, • ■ ■ , L2 n/?b0 J - 1} := 5 b0 
and generate a partition of tS b0 randomly by independently assigning every index Wbo G 5 b o to a 
set iS b o,i, with a uniform distribution over the indices i G {0, . . . , [2 nifbl J — 1} := Sbi- We denote 
by Sbo(wbo) the index i of <S b o,i to which u; b o belongs. 

3) Phase 3: Generate random (n • A3 n )-length sequences 

. u ( a 3 \w a2 ) i.i.d with p^(u a ), w a2 G {0, 1, • • • , L2 nif -J - 1} := S a2 
4 This choice of A3 , a a , oib is on the boundary of the closure of the achievable rate region. 



y r (1) (Wao) i-i-d. With p«(£ r ) = E^^pW^rbr) , ^a0 G {0, 1, • • ■ , [2 nR >» \ - 1} := S a 
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. u^(w b2 ) i.i.d with pW(u b ), w b2 G {0, 1, ■ ■ ■ , [2 nR ^\ - 1} := S h2 
and define bin Bj := {w a2 \w a2 G [(j - 1) ■ \2< R >*- R ^\ + 1, j • L2 n ( Ra2 -- Ral )j]} for j G 5 ai . 
Likewise, C k := {w b2 \w b2 G [(k - 1) • |_2 n ( J * b2 ~ ilbl )j + 1, jfe . L2 n ( ft > 2 - i? «)j]} for fc G S b i- 
4) Phase 4: Generate random (n • A4 n )-length sequences 

. x[ 4) (w a0 ,u;bo) i-i-d with p( 4 )(rr r ), w a0 G S a0 and u; b0 G <S b0 . 
Encoding: During phase 1 (resp. phase 2), the encoder of node a (resp. b) sends the codeword x a (w a ) 
(resp. x. h 2 \w b )). At the end of phase 1, relay r compresses the received signal y? 1 ' into the message 
w 3 q if there exists a w a o such that (y£ , yr (^ao)) G A^(Y r Y r ). Similarly, r compresses yi into the 
message w b Q at the end of phase 2. There exist such w 3 q and io b o with high probability if 

R a0 = Ai, n /(Y r (1) ; y r (1) ) + e (66) 
i? b0 = A 2i?t /(y r (2) ;y r (2) ) + e (67) 

and n is sufficiently large. Also we choose R a \, R b i, R a2 and R b2 as: 

i? a i = a a R a0 = a a (A lirt /(y r (1) ; y r (1) ) + e) (68) 
Rbi = n b R b0 = a b (A 2 , n I(y r (2) ;y r (2) ) + e) (69) 

and 

i?al < R 3 2 = A 3 ,n/(^i 3) \ Y® ) - 4e (70) 

Rbi < Rb2 = A 3 ,„/(^ 3) ;y a (3) ) - 4e. (71) 

From the code constructions of w a2 and w b2 , R a ± and i? bl have to be less than R a2 and R b2 , respectively. 
Then the relay constructs w a \ = s a o(w a o) and w b \ = s b o(w b o). To choose w a2 and w b2 , the relay first 
selects the bins B Wal and C Wbl and then it searches for a pair (w a2 ,w b2 ) G B Wal x C Wbl such that 
(u a 3 \w a2 ), u b 3 \w b2 )) G A^(U a U b ). Such a (io a 2,io b 2) exists with high probability if 

R al + R bl < R a2 + R b2 - A 3 ,„/(^ a (3) ; ) - e' (72) 

from the Lemma in [6]. The relay then sends generated i.i.d. according to p( 3 \x r \u a , u b ) with 
u^\w a2 ) and u h 3 \w b2 ) during phase 3. Finally, the relay sends x^(to a o, u; b o) during phase 4. 

Decoding: Node a decodes w b2 after phase 3 using jointly typical decoding. Then a estimates w b \ 
from the bin index of w b2 . Node a decodes w b Q if there exists a unique -oj b o such that w b Q G <Sb0i5 b i> 
(x[ 4) (u; a o,* b o), yi 4) ) G AW(X r Y a ), (x a 1) ( U ; a ) ) y«(^ao)) G A^{X a Y r ) and (y[ 2) («; b0 ), y^ 2) ) G A^(Y r Y a [ 
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(21 (2) (2) 

. After decoding w^o, node a decodes w b using jointly typical decoding of the sequence (x b , y r (w b o) > Ya )• 
Similarly, node b decodes u) a . 
Error analysis: 

P[E b , a ] < P[E$ U E$ U P^] (73) 
< P[E$] + P[P r (4) |P r (3) ] + P[E$ \E® n (74) 

We define error events in each phase as follows: 

1) E$$ = Eif 1 U Eif 2 U E^ 3 U P r (3) ' 4 U E^f 5 . 

E^' 1 : there does not exist a u; a0 such that (y r (1) ,y r (1 W)) e A^(Y r Y r ). 
Er 3 ]' 2 : there does not exist a u; b0 such that (yP'.yPVbo)) G ^(yy). 

-Er,a^' 3 : there does not exist a pair (w a2 , w b2 ) G £ Sa0 (™ >0 ) xC Sbo(lUb0 ) such that (u a 3) (iw a2 ), u b 3) (w b2 )) G 

^ (3) (^af/b). 

4fa } ' 4 : (u£\w b2 ),y®)?AW(U b Y a ). 

P r (2) ' 5 : there exists w b2 + w b2 such that (u b 3) (u) b2 ), yf ] ) G A( 3 )([/ b y a ). 

2) 4 4 a } = 23gP U Eif 2 U 4 4 a } ' 3 U 23<2' 4 . 

P^' 1 : there does not exist a u; a0 such that (y r (1) ,yr (1 W)) G A«(yy). 



^r,a 

p(4),3 
^r,a 

^r,a 



there does not exist a u? b0 such that (y r (2) ,y r (2) (w b0 )) G A( 2 >(y r £). 

(x[ 4) Ko,^bo),ya 4) )^^ (4) (x r y a ). 

there exists (w a0 , w b0 ) where w a0 ^ w a0 and w b0 / iu b0 such that (x[ 4) (w a0 , w b0 ), y a 4) ) G 
A( 4 )(X r y a ), (x^K),^ 1 ^)) G AW(X a y r ), (y r (2) (^bo),yi 2) ) G A®{Y t Y a ), and 



*b0 G S bQiSbo ( Wbo y 

E^' 5 : there exists w b0 ^ w b0 such that (x[ 4) (w a0 , w b0 ), y^ ) G A< 4 ) (X r Y a ), (y^ 2) (w b0 ), y (2) 
A^(Y r Y a ), and * b0 G S b0 , Sbo{Wbo) 

3) ^=4ru^. 

n (4),l . / (2), v (2) A (2). 



a ;G 



<T : (x b 2j K),yi 2J ,y^Ko)) 0A( 2 )(x b y a y r) 

pg' 2 : there exists w b + w b such that (x { b 2) (w b ) , y a 2) , y^ (w b0 )) G A^(X b Y a Y r ). 

Then, 



P[E$] <P[E^ A ] + P[P r (3) ' 2 ] + P[E®> 3 ] + P[P r (3) ' 4 ] + P[£r ( ?a' 5 ] (75) 

<3e + P[,D( 3 )(u b (u; b2 ),y a )] +P[U lib2 ^ b2 Z)( 3 )(u b (TZ; b2 ),y a )] (76) 
<4 e + 2 ™(^2-A 3 .„/(C/ b (3, ;y a (3) )+3e) (77) 
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In (/J), P[£$ ),:L ], PfPrg' 2 ] and P[E$' 3 ] are less than e due to <[66]>, (HI]) and (72]), respectively. 

P^Sl^g?] ^P^S' 1 !^] + P[E$$*\E®] + P[E$*\&J]+ 

P[M2*\M$]+P[B$$ fi \Mi] (78) 
<P[5W(x r ( Wa0 ,w b0 ),y a )]+ 



P 



u a a0 #» a o P (4) (x r (u) a0 , w b o ) , y a ) , P (1) (x a (w a ) , y r (u) a0 ) ) , 



P (2) (y r < 



Wb0),ya),Sb0(^b0) = ™bl 



+ 



P[Um bo ^ bo r^ (4) (x r (w a o,Wbo),ya),^ (2 Hyr(*bo),ya),Sbo(*bo) = Wbl] (79) 

<£ + 2 n(fi, +i? b o-A 4 .,a(X< 4 ^y (4> )-A 1 ,^ 

2 n( J R b o-A 4 ,„/(X r (4) ;y a (4) )-A 2 ,„/(y r (2) ;y r <2> )-a bJ R bo +e'") ^ 

In C§), P[E^} ,l \E^}] and Pfpg^lP^J] are zero since pgP and P r (4) ' 2 are the same as P®' 1 and 
Pg , respectively. In (f80b . the bound for Pbo in the second term is implied by that in the third term 
since P a0 - A lin /(y r {1) ;X a {1) ) = A lin /(Y r (1) ; Y r (1) |X a (1) ) + e > 0. 

P[*$\s®nM®] <p[eIT\e% npg } ] + P[Pi 4 a ),2 |pg ) npgi] (81) 

=P[D( 2 '(x b ( Wb ),y b ,y r (» b0 ))] + P[U^ b ^ b P( 2 ) (x b (to b ), y b , y r («>bo))] (82) 

< e + 2 n(%-A 2 ,„7(x b <2) ;y r <2) ,y <2) )+3 e ) (83) 

Since e > is arbitrary, a proper choice of a b , the conditions of Theorem [l2j (T7TT >. and the AEP 
property guarantee that the right hand sides of (1771 ). (l80l ) and (l83l ) corresponding to (1711 ). (l60l ) and (l55T ) 
vanish as n — > oo. Similarly, P[P ajb ] — > as n — ► oo. By the Caratheodory theorem in [9], it is sufficient 
to restrict \Q\ < 13. ■ 

When the h a and /i b links are of different strength, a scheme in which one link uses CF and the other 
uses DF may provide a larger rate region than if both links use CF. In the next theorem, we provide a 
rate region for a TDBC scenario in which the forward link uses DF and the reverse link uses CF. 

Theorem 14: An achievable rate region for the half-duplex bi-directional relay channel with a mixed 
TDBC protocol, where the a — ► r — > b link uses decode and forward and the b — > r — > a link uses 
compress and forward, is the closure of the set of all points (P a , P b ) satisfying 

R a <min{A 1 I(X a (1) ; y r (1) ),A 1 I^^ (84) 

P b < A 2 I(xf ;# 2) ,Y a (2) |Q) (85) 
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subject to 

A 2 /(y r (2) ; y r (2) |ya (2) ,Q) < mm{A 3 I(U®,U®;Y®\Q), A 3 /(C/f ) ;C/ r (3) ,n (3) |Q)} (86) 
over all joint distributions, 

p(q, x a ,x b , x n u b , u r , y a ,y b ,y r , y r ) = p{q)p W {x a ,y b ,y r )p^ (x h , y a , y r , y r \q)p {3 \u b , u r , x r , y a ,y b \q) 

(87) 

where 

P {l) {x a ,y h ,y r ) =p {1) (x a ) P ^(y b ,y r \x a ) (88) 

P {2 \x bl y a ,y r ,y r \q) = p {2) (x b \q)p {2 \y a ,y r \x b )p ( - 2 \y r \y n q) (89) 

P i3 Hu a , u b , u r , x n y a ,y b \q) = p {3 \u b , u r \q)p^ {x r \u b , u r , q)p {3) {y a ,y b \x r ) (90) 

with | Q\ < 6 over the alphabet X a x X b x X x x U b x U x x >^ a 2 x ^ x ^ r 2 x j) r . □ 
Remark 15: We use random binning and a Gel'fand-Pinsker coding scheme in Theorem [T4j The 
detailed proof is provided in Appendix U 

In contrast to the MABC protocols, in the TDBC protocols, the mixed TDBC protocol is not outer 
bounded by the DF TDBC protocol. In the TDBC protocol, each terminal node obtains side infor- 
mation, used during decoding, when the opposite node transmits (phase 1 for b and phase 2 for a). 
The data rate as well as the phase durations Ai and A 2 of the TDBC protocol may vary, and we 
note that one can easily find cases in which the mixed TDBC protocol outperforms the DF TDBC 
protocol. For example, suppose the channel of the link r <-» a is good enough such that I{Xr>-Y a yo >) > 
mm{I(U^ 3 \u^ 3) ;Y a {3) )J(U^ ) ;U^\Y a {3) )} » and the link mutual informations between b <-> r and 
b <-> a are the same such that I(X {2) ]Y r {2) ) = I(X {2) ;Y a {2) ). Furthermore, take A a = e lt A 3 = e 2 , 
where €\ and e 2 are small positive numbers. We assume the input distributions are fixed, i.e. \Q\ = 1. 
Then, in the DF TDBC protocol, from Theorem |4j an achievable rate region is given by: 

R a < e' (91) 

#b < (1 - ei - e 2 )I{X {2) -Y a {2) ) ~ I(xi 2) ;Y a {2) ), (92) 

for some e'. In the mixed TDBC protocol, we find a choice for min{/(t/ r (3) , ; Y a {3) ),I(U^ 3) ; f/ r (3) , Y a {3) ) } 
which satisfies (l86l ) for the same Aj's (this is possible as the r «-> a channel is very strong). Then we 
obtain the following achievable rate region from Theorem [14] : 

R a < e" (93) 

R b < (i - ei - 6 2 )/(xf ; y a (2) ,y (2) ) « /(xf ; y a (2) ) + /(x£ 2) ; y r {2) |y a {2) ), (94) 
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for some e". By properly choosing ei and the regions bounded by d9~TT ), d92l ), (|93l ) and (|94l ) demonstrate 
the case in which a larger rate region is achieved through the mixed TDBC protocol. 

IV. Gaussian Case 

We now assume all links in the bi-directional relay channel are subject to independent, identically 
distributed white Gaussian noise. The commonly considered Gaussian channel will allow us to visually 
compare different achievable rate regions for the bi-directional relaying channel. Definitions of codes, rate, 
and achievability in the memoryless Gaussian channels are analogous to those of the discrete memoryless 
channels. 

We apply the previous results to the Gaussian channel. Since strong typicality does not apply to 
continuous random variables, the achievable rate regions from the theorems in the previous section do 
not directly apply to continuous domains. However, for the Gaussian input distributions and additive 
Gaussian noise which we will assume in the following, the Markov lemma of [17], which generalizes 
the Markov lemma to the continuous domains, ensures that the achievable rate regions in the previous 
section hold in the Gaussian case. We use Gaussian input distributions since we assume the average 
power constraint. 

The corresponding Gaussian channel model is: 

Y 3 [m] = h r3 X r [m] + h b3 X b [m] + Z 3 [m] (95) 
Y b [m] = h rb X r [m] + h ab X a [m] + Z b [m] (96) 
Y r [m] = h 3r X 3 [m] + h br X b [m] + Z r [m] (97) 

where X a [m], X b [m] and X r [m] follow the input distributions X { 3 £) ~ CW(0, P a ), X^ ~ CM{0,P b ) 
and X^ ~ CM(0, P r ) respectively during transmitting, where m G [n X^j=o A?> + ^ n S^=o A?>] anc ^ 
CJ\f(fi, a 2 ) denotes a complex Gaussian random variable with mean n and variance a 2 , and t corresponds 
to the appropriate phase. If node i is in transmitting mode, the transmit power is bounded by Pi, i.e., 
^[^ 2 ] < Pi- If node i is in receiving mode, the input symbol does not exist in the above mathematical 
channel model. For example, in the first phase of the TDBC protocol, the corresponding channel model 



is : 

Y b [m] = h 3b X 3 [m] + Z b [m] (98) 

Y r [m] = h„X a [m] + Z r [m]. (99) 
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TABLE IV 
Input and output distributions 



Dist. of 


during TX 


during RX 


x^ 

7 W _ 7 W _ 7 W 

Zi 3 — Z/ b — Zj, 

yW 

Z^ 


CAf{0,P 3 ) 
CAf(0, P b ) 
CA/"(0, P r ) 

N/A 

N/A 

N/A 


N/A 
N/A 
N/A 
CA/"(0, 1) 

c^(o,pf) 
CAr(o,pf 



In the above fo^ is the effective channel gain between transmitter % and receiver j, which is modeled as 
a complex number. We assume that the channel is reciprocal such that hij = hji and each node is fully 
aware of h ar , h^ r and h a b (i.e. we have full CSI). The noise at all receivers Z a , Z^, Z r is of unit power, 
additive, white Gaussian, complex and circularly symmetric. For convenience of analysis, we also define 
the function C(x) := log 2 (l + x). 

* (£) 

For the analysis of the Compress and Forward scheme, we assume Y T are zero mean Gaussians and 
define P<f } := E[(Y {£) ) 2 ] , pf> := E{(Y (l) ) 2 } and a$ := E[Y {e) Y {e) ]. Then the relation between the 
received Y r [m] and the compressed Y r [m] are given by the following equivalent channel model: 

Y[m] = h ri [m]Y[m] + Z f [m] (100) 

where Y r [m], Y r [m] and Z f [m] follow the distributions F r W ~ CAf(0, Pf ] ), Y r W ~ CN{<d,Pf ] ) and 
zf ~ CN^Pf - i^f) and /i rf [m] = ^J, where m € [n£jlj A j,n + l,n^ =0 A,- n ]. We note 
that in the following, P~ and Oj) are unknown variables corresponding to the quantization which we 
numerically optimize. 

We consider four different relaying schemes (i.e. ways in which the relay processes and forwards the 
received signal) for each MABC and TDBC bi-directional protocol: Amplify and Forward (AF), Decode 
and Forward (DF), Compress and Forward (CF), and Mixed Forward (Mixed). In addition to achievable 
rate regions, we apply outer bounds of the MABC and TDBC protocols to the Gaussian channel. 

A. Amplify and Forward 

In the amplify and forward scheme, all phase durations are equal, since relaying is performed on a 
symbol by symbol basis. Therefore, Ai = A 2 = 5 for the MABC protocol and Ai = A 2 = A 3 = | 
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for the TDBD protocol. Furthermore, relay r scales the received symbol y r by J -§- to meet the transmit 

V y 

power constraint of P r . The following are achievable rate regions for the amplify and forward relaying: 
. MABC Protocol 

l c ( \h 3t \ 2 \h hr \ 2 P a P r 

2 Vl^ar| 2 Pa + |^br| 2 ^ + |^br| 2 Pr + l 



R > < ,C [ n - | -;, 2 , ) (101) 



t> A C ( \h 3r \ 2 \h br \ 2 P b P r 

Rb< 2° I \h^P a + \h br \ 2 P b + \h ar \*P r + l ] (1 ° 2) 



. TDBC Protocol 



^<l c h 2p - + MirMMM^ ) <io3) 

B. Decode and Forward 

Applying Theorems [3] and 0] to the Gaussian case, we obtain the following achievable rate regions: 
. MABC Protocol 

R 3 < mm{A 1 C(\h 3r \ 2 P 3 ),A 2 C(\h br \ 2 P r )} (105) 

R b < min{AiC(|/i br | 2 PL), A 2 C{\h 3r \ 2 P r )} (106) 

R 3 + R b < A!C(|/i ar | 2 P a + \h br \ 2 P b ) (107) 

. TDBC Protocol 

P a < min{AiC(|h ar | 2 P a ), A ± C \\h 3b \ 2 P a ) + A 3 C(|/i br | 2 P r )} (108) 

R b < mm{A 2 C(\h br \ 2 P b ),A 2 C(\h 3b \ 2 P b ) + A 3 C(|/i ar | 2 P r )} (109) 

When obtaining the regions numerically, we optimize A/s for the given channel mutual informations to 
maximize the achievable rate regions. 

C. Compress and Forward 

Applying Theorem 151 and [121 to the Gaussian case, we obtain the following achievable rate regions: 
. MABC Protocol 

P a <A 1 cf a K ))2|/l m 2jPa m I ( 11Q ) 

* < AlC I (^MP, 

I P{, {Py ) 2 ~ ( a i ) 2 (Py ~ 1) 
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where, 



mm < 



C{\h br \ 2 P r ) 



c 



{<T™n\h,\*p,+i) 
p?Hppy-{oPyp?> 



+ C{\h hr \ 2 P r ) c 



= \h a r\ 2 P a + \h br \ 2 P b + l 



C(\h 



, 2 P) 

ar | 1 r) 



(4 1) ) 2 (|fe br | 2 P b +l) 



+ C(\h ar \ 2 P r ) 
(112) 
(113) 



TDBC Protocol 

One can show that (l56l ) - (I58T ) Marton's bound is equivalent to the capacity region of the Gaussian 

broadcast channel with Costa's setup as follows: let \h ra \ > \h rb \ and we set 

t 

U b \m\ = V r [m] + a?7 a [m] 

Y a [m] = h ra (V r [m] + U 3 [m\) + Z a [m] (114) 
_ nH = h rh {V r [m] + U a [m}) + Z b [m] 

where V r [m] and U a [m] follow the distributions V r (3) ~ CM(0,(3P r ), £/ a (3) ~ CJ\f(0, (1 - /3)P r ) 
respectively during phase 3, m 6 [n(Ai n + A2, n ) + 1, i], where (0 < j3 < 1) and ^[V - / 3 ^?/^] = 0, 
(3) rr( 3 ) ««. „w ii m ™. ^ „ - 



In phase 3 < 



i.e., V r , C^a are independent. Also we take a — , h \2p P+ i 

r(3). v( 3 ) 



. Then 



r , |Ah,1 3 (1-/?)pA 

1 |/l,b| 2 /9Pr+l J 



(115) 



-/(t/a (3) ;^ 3) ) = C(|^ ra | 2 /3P r ) 
Similarly, we obtain the bounds in the case \h ra \ < \h rb \. These are the same as the capacity region 



of the Gaussian broadcast channel ((15.11) and (15.12) in [4]). 



R a < AiC \h ab \ 2 P a + 



R h < A 2 C \h ab \ 2 P b + 



(2)x 2 | 



(a"y'r\h br \ 2 P b 



pf(pi 2) ) 2 -(4 2) ) 2 (^ (2) 



1) 



(116) 



(117) 
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where, 

/ (a {1) ) 2 \ 

If \h n \ < \h rb \ : a^C [ ^ m J m y L flKo I < A 3 C (/J^rbl'Pr) (118) 

I <A 3 cf (1 "^ ra|2Pr V (H9) 




,3|/l ra | 2 P r + l 

( ( !„7?o^ rb|2 f r ) 



/?|/l rb | 2 Pr + l 

a b A 2 C I ^ 7 I < A 3 C (/?|/* ra | 2 P r ) , (121) 



and 



p^p^-^yj \pf\p^?-{^?p^ 

br | 2 P r ) + AiC ( (4 1) ) 2 |^ab| 2 |/lr,| 2 P, 



< A 3 C(|h ar | 2 P r ) + A 2 C f _ (4 2) ) 2 l^b| 2 |^b| 2 Pb 

\(Pj } ) 2 Pf (l^abPPb + 1) " (4 2) ) 2 ^ab| 2 |^rb| 2 Pb 

P y (1) = |/iar| 2 Pa + l (124) 

Pf) = |/ibr| 2 Pb + l (125) 

< a a ,a b ,/3 < 1 (126) 

Again, when numerically obtaining the regions, we optimize P~ , a y , Ai, a 3 , a b and p to maximize 
the region boundary. 

D. Mixed Forward 

Applying Theorem [8] to the Gaussian case with Costa's setup in [3] we have the channel in the relay 
broadcasting phase in the Mixed MABC protocol as: 

U r [m] = V r [m] + aU b [m] 
In phase 2 <( Y a [m] = h ra (V r [m] + U b [m]) + Z a [m] (127) 
YbH = h rb (V r [m] + U b [m\) + Z b [m] 
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where V T [m] and U b [m] follow the distributions V r {2) ~ CM{0,(3P r ), U { h 2) ~ £/V(0, (1 - /3)P r ) during 
phase 2, m G [Ai, n -n+l,n], where (0 < (3 < 1), and P[K (2) f7 b (2) ] = 0, i.e., V r {2) , U {2) are independent. 

Similarly, we construct the channel in the relay broadcasting phase in the Mixed TDBC protocol. Then 

(£) (£) 

we obtain the following achievable rate regions, where we numerically optimize a, [5, P~ , &y and A^ 
to maximize their boundary. 
. MABC Protocol 



\h 3 r\ 2 P 3 \ a . ( /?P r (|/l rb | 2 P r + l) 

R 3 < mm <^ AiC - , A 2 log 1 



:)} 



|/tbr| 2 Pb + I/ ' 62 Vl^b| 2 (l-a) 2 /3(l-/3)P r 2 + /3P r + a 2 (l-/3)P r 

(128) 

m^ ))2 y ^ m I 

where, 

(^^(l^pPb + l) \ < ( A2C , ( | hra |2p r)) ( |^ a |2 (1 _ a)2(1 _ ^ + « 2 (1 



(130) 

P y W = |/iar! 2 Pa + |/jbr| 2 Pb + l (131) 

. TDBC Protocol 

P a < min{AiC(|/i ar | 2 P a ), 



AiC(|/ lab | 2 P a ) + A 3 log 2 



(3P r (\K b \ 2 P r + 1) 



|/t rb | 2 (l - a) 2 /3(l - /3)P 2 + /?P r + a 2 (l - (3)P t , 

(132) 



Pb < A 2 C I |/l a b| 2 Pb H TTvj 77yT~ V - — TZ\ 77^ — I (133) 

pf (pf) 2 -(4 2) ) 2 ^ (2) + (4 2) ) 2 



where, 



P 9 ? 1 <min(A 3 C(|/ lra | 2 P r ),A 3 cf|/ lra | 2 (l-a) 2 (l-/3)P r + a2(1 # 



^i 2) ^ 2) -(4 2) ) 2 / l v 1 v ; v ' 

(134) 

= |/l b r| 2 Pb + 1 (135) 
n * |^ab| 2 Pb |^rb| 2 Pb 



|/lab| 2 Pb + l |/l rb | 2 Pb + l 



(136) 
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E. Outer Bound 

Applying Theorems Q] and [2] to the Gaussian case, we obtain the following outer bounds. We optimize 
A/s for given channel mutual informations to maximize these outer bounds. 
. MABC Protocol 

R a < minlAxCd^pPa), A 2 C(|/ lbr | 2 P r )} (137) 
R b < mm{A 1 C(\h br \ 2 P b ),A 2 C(\h ar \ 2 P r )} (138) 

. TDBC Protocol 

R a < mmiAMKrfPa + \h ab \ 2 P a ), &iC(\h ab \ 2 P a ) + A 3 C(|fc br | 2 P r )} (139) 
R b < min{A 2 C(|/i br | 2 P b + \h ab \ 2 P b ), A 2 C (\h ab \ 2 P b ) + A 3 C(|/i ar | 2 P r )} (140) 
R a + R b < AiC(|/i ar | 2 P a ) + A 2 C(|/i br | 2 P b ) (141) 

V. ACHIEVABLE RATE REGIONS IN THE GAUSSIAN CHANNEL 

In order to obtain an intuitive feel for the regions and to illustrate that the regions are not subsets 
of one another, the bounds described in Section JV] are plotted in this section for a number of different 
channel configurations. We first compare the rate regions obtained by the bi-directional protocols and 
outer bounds in cases in which the links are symmetric (h ar = h br = 1, h ab = 0.2) as well as asymmetric 
(h ar = 0.6, h br = 20, h ab = 0.5 and h ar = 20, h br = 0.6, h ab = 0.5) for two different transmit SNRs of 
and 20dB. The protocols considered are: 

1) AF MABC: Amplify and Forward Multiple Access Broadcast. 

2) AF TDBC: Amplify and Forward Time Division Broadcast. 

3) DF MABC: Decode and Forward Multiple Access Broadcast. 

4) DF TDBC: Decode and Forward Time Division Broadcast. 

5) CF MABC: Compress and Forward Multiple Access Broadcast. 

6) CF TDBC: Compress and Forward Time Division Broadcast. 

7) Mixed MABC: Mixed scheme with Multiple Access Broadcast. 

8) Mixed TDBC: Mixed scheme with Time Division Broadcast. 

9) Outer MABC: Outer bound when using MABC protocol. 
10) Outer TDBC: Outer bound when using TDBC protocol. 

We then proceed to examine the maximal sum-rate R a + R b of the ten schemes as a function of the 
transmit SNR. Finally, we evaluate the maximal sum-rate and maximal constrained sum-rates (that is we 
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require the rate R a = R b as well as R a = 2R b ) of the schemes as a function of the relay position. The 
main conclusions to be drawn are that different schemes are optimal under different channel conditions. 
We provide further discussions in the following subsections. 

A. Achievable rate region comparisons 

We compare the achievable rate regions and outer bounds of the 10 aforementioned protocols for both 
symmetric and asymmetric source to relay channel gains at transmit SNRs of and 20dB. 

1) Symmetric Case: In this case h ar = = 1 (Figs. E). In the low SNR regime, the DF MABC 
protocol dominates the other protocols. The MABC protocol in general outperforms the TDBC protocol 
as the benefits of side information and reduced interference are relatively small in this regime. The DF 
scheme outperforms the other schemes since the relatively large amount of noise in the first phase (and 
the second phase in the TDBC protocol) can be eliminated in the DF scheme, which cannot be done 
using the other schemes. In contrast, the DF TDBC protocol dominates the other protocols at high SNR 
since the direct link is strong enough to convey information in this regime. 

In the high SNR regime (when P a = P b = P r = P is sufficiently large), the AF MABC protocol 
outperforms the DF MABC protocol. From (11011) . (11021) . the achievable rate region of the AF MABC 
protocol is: 

If P 2 \ 1 
R a < - log 1 + ~ - log P (142) 

If P 2 \ 1 
flb<- loglH w-logP (143) 

also from dl05l) , dl06b and d 107b . the achievable rate region of the DF MABC protocol is: 

R a <min{Ai,l- Ai}-log(l + P) « min{Ai, 1 - Ai} • log P < ^logP (144) 

P b <min{Ai,l- Ai}-log(l + P) « min{Ai, 1 - Ai} ■ log P < ^logP (145) 

R a +R b < Ai log(l + 2P) « Ai log P (146) 

From (11461 ). we can conclude that the achievable rate region of the DF MABC protocol is outer-bounded 
by the AF MABC protocol. 

In the TDBC protocol, the CF scheme does not outperform the DF scheme since the DF uses two 
parallel channels in phase one and three, while the CF uses one channel in phase one with two receivers. 
In other words, Pf F < AiC(-) + A 3 C(-) for the DF as opposed to R° F < AiC(£ •) for the CF 
scheme. However, under the MABC protocol, the CF scheme outperforms the DF scheme in the high 
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SNR regime. This is because the interference of the transmission of two terminal nodes affects the DF 
MABC scheme due to the multiple-access nature but not the CF scheme (as it must not decode the 
signals). 

The achievable rate region of the Mixed TDBC protocol lies between the CF TDBC protocol and the 
DF TDBC protocol. In the TDBC protocol, max Rb R^ IX = max Rb R® F , where Rf IX is the data rate 
of node a in the Mixed scheme. Here the max is taken over all rates in the achievable rate regions. The 
max R a IIX is achieved by taking A2 = 0. The rate R® F is similarly defined and max R^ F is achieved 
in an analogous manner. Therefore, the point (maxfif F , 0) lies in both the Mixed scheme and the DF 
scheme. While in the TDBC scheme a particular rate may be set to by indirectly setting the appropriate 
interval A, to 0, in the MABC protocol this is not possible. In the MABC protocol, even when R 3 = 
or i?b = 0, the transmit power (P a or P^) remains constant in the first phase and does not decrease to 
0, acting as the additional noise for the opposite transmission. Therefore, max^ b 

R MIX < 

max/j b i?a . 

This interference seen in MABC protocols is especially pronounced in the high SNR regime, where 
the gap between the intercept points of the Mixed scheme and the DF scheme is seen to grow as the 
interference increases (with SNR). We note that this effect is due to our assumption that both transmitters 
transmit at full power regardless of the rate. If we were to allow for optimization of the transmission 
power, larger achievable rate regions for the Mixed MABC protocol could result. 

In Figs. [6] and [71 the AF scheme is always outer bounded by the CF scheme. We thus expect that when 
relay r does not know the codebooks of a and b (and hence cannot decode as in the DF scheme), that 
the CF scheme is a better strategy than the AF scheme. 

In the low SNR regime, the achievable rate region of the DF MABC protocol and the outer bound 
of the MABC protocol are tight, while in the high SNR regime, the achievable rate region of the CF 
MABC protocol is tight. For the TDBC protocol, there is a very small gap between the achievable rate 
region of the DF TDBC protocol and the outer bound of the TDBC protocol since interference is not an 
issue for the TDBC protocol and hence decoding is thus, intuitively, near optimal. 

2) Asymmetric Cases: In these cases h ar = 0.6, h^ x = 20, h a ^ = 0.5 (Figs. [8j [9]) and h ar = 20, h^ r = 
0.6, /i a b = 0.5 (Figs. ITOl ITTb . Note that these two asymmetric cases are different for the mixed forwarding 
cases, which assume CF in one direction and DF in the other. In the low SNR regime, the CF TDBC 
and mixed TDBC protocol achieve the best performance in Fig. [8] and Fig. [TOj respectively. However, in 
the high SNR regime, the DF MABC protocol and the DF TDBC protocol yields larger regions than the 
other protocols. In contrast to the symmetric case, the AF MABC protocol is not outer bounded by the 
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Fig. 6. Comparison of bi-directional regions with h 3r = Fig. 7. Comparison of bi-directional regions with h a , = 
h b , = 1, h 3b = 0.2, P a = P b = P = dB and iV a = iV b = h b , = 1, h sb = 0.2, P a = P b = P r = 20 dB and N a = 
N r = 1. N b — N r = 1. 



CF MABC protocol. 

The mixed forwarding scheme is the only one which has different performance in the two asymmetric 
cases. In the mixed MABC protocol, if h^ r > h ar then the noise seen at relay r when it decodes w a is 
larger than when h\, r < h ar . Therefore, the corresponding achievable rate region is also relatively smaller. 
In particular, in the high SNR regime (Fig. [9]), the achievable rate region for the mixed MABC protocol 
is only able to achieve, under our input assumptions, rate R a close to because of the effective noise 
from node b during the first phase. In the mixed TDBC protocol, if h ar > h^, then we have a larger 
achievable rate region since the first link is more critical to the performance of the DF scheme. As the 
SNR increases, the difference between the two asymmetric cases decreases. 

B. Maximum Sum Data Rate 

In this subsection we plot the maximum sum-rate R a + i?b as a function of the transmit SNR for 
the symmetric and two asymmetric cases of the previous subsection. As expected, different schemes 
dominate for different SNR values. The sum-rate is basically proportional to the SNR in dB scale since 
the sum-rate is roughly the logarithm of the SNR. In Fig. [12] around 12 dB the relative performance of 
the CF MABC protocol and the DF MABC protocol changes. At lower SNRs, the DF MABC protocol 
is better, while at higher SNRs, the CF MABC protocol is better. The AF MABC protocol is always 
worse than the CF MABC protocol in the symmetric case (Fig. W2\ . In the TDBC protocol, the sum-rate 
of the mixed TDBC protocol lies between the DF scheme and the CF scheme in Fig. [12] 
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C. Relay position 

In this subsection we plot the maximum sum-rate R 3 +Rb as a function of the relay position d ar = (d a b 
(0 < C < 1) when the relay r is located on the line between a and b. Thus, d^ = (1 — C)d a b- We apply 
/i a b = 0.2 and P a = P^ = P r = 20 dB and let \hij\ 2 = k/dfj 8 for k constant and a path-loss exponent 
of 3.8. We consider three constraints on the sum-rate in the three Figs. [T5J [16] and [T7] In the first, the 
sum-rate is maximized without any additional constraints. For the latter two we consider more realistic 
scenarios in which the sum-rate is constrained. In many communication systems, uplink and downlink 
rates are not equal. More specifically, it is not uncommon for the dowlink rate to be 2 to 4 times greater 
than that of the uplink. In Figs. [16] and [T7] we plot the maximal sum-rate of the protocols under a 
o = i? a /i?b rate ratio restriction. These constrained sum-rates are obtained by optimizing the Gaussian 
regions of Section [TV] with the additional constraint R a = aR^. 

For the MABC protocol, if the relay location is biased (not midway between nodes a and b), then 
the DF MABC protocol outperforms the CF MABC protocol and the AF MABC protocol outperforms 
the CF MABC protocol. This effect is more explicit in the constrained cases. In contrast, for the TDBC 
protocol, in order of increasing complexity, (and performance), the relaying schemes are AF, CF, Mixed 
and DF in all cases when (0.1 < ( < 0.9). As expected, the sum-rate for the mixed MABC protocol is 
worse than those of the other protocols. 

The sum-rate plot for the mixed protocol is not symmetric since it uses different forwarding schemes 
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Fig. 10. Comparison of bi-directional regions with h a , = 
20, h br = 0.6, /lab = 0.5, P a = P b = P = dB and 

iVa = iVb = N, = 1. 




R 

a 



Fig. 11. Comparison of bi-directional regions with h 3r = 
20, h br = 0.6, h 3b = 0.5, P a = Pb = Pr = 20 dB and 
N 3 = N b = N r = 1. 



for each link. In addition, the sum-rates in the constrained case where a = 2 (Fig. ITTb are asymmetric even 
for non-mixed protocols. The intuitive reason for this is that the rates are constrained in an asymmetric 
way, and hence a particular, non-midpoint distance will be optimal even for CF and DF forwarding 
schemes. The performance of the a = 2 sum-rate for the DF MABC protocol is remarkably asymmetric 
where it peaks and almost touches the outer bound at Q = 0.6. These plots and region optimizations may 
be useful when determining the optimal relay position subject to particular rate constraints. 

VI. Conclusion 

In this paper, we derived achievable rate regions for 4 new half-duplex bi-directional relaying protocols. 
We have specialized the 8 different achievable rate regions and the 2 outer bounds to the Gaussian case 
and numerically evaluated them under various channel conditions. For the MABC protocol, DF or CF 
is the optimal scheme, depending on the given channel and SNR regime. In the TDBC protocol, the 
relative performance of the forwarding schemes depends on the given channel condition. Notably, we 
have determined an example of a channel condition in which the mixed TDBC protocol outperforms the 
other proposed protocols. In general, the MABC protocol outperforms the TDBC protocol in the low 
SNR regime, while the reverse is true in the high SNR regime. 
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Fig. 12. Maximum sum-rate of the 8 bi-diretional protocols and 2 outer bounds at different SNR. Here h ar — hbr — 1 and 
/lab = 0.2 



Appendix I 
Proof of Theorem [141 

Proof: Random code generation: For simplicity of exposition, we take \Q\ = 1. 

1) Phase 1: Generate random (n • Ai >n )-length sequences 

. x^K) i.i.d. with pM(x a ), w a G S a = {0, 1, • ■ ■ , [2 nR -] - 1} 
and generate a partition of S a randomly by independently assigning every index w a G S a to a set 
<5 a ,i, with a uniform distribution over the indices i G {0, . . . , L2 nifa0 J — 1} := S a o. We denote by 
•Sa(^a) the index i of S a ^ to which ui a belongs. 

2) Phase 2: Generate random (n • A2, n )-lengfh sequences 

. x[, 2) (u/ b ) i.i.d. with ^(^b), %G5 b = {0,l,---, L2 ni?b J - 1} 

. yP } Ko) i.i.d. with P W(y r ) = J2y r p {2) (yr)p {2) (y r \yr) , w r0 e {o, l, • ■ ■ , [2 nR ->\ - 1} : = 5 r0 

3) Phase 3: Generate random (n • A3 n )-length sequences 

. u[ 3) (-u; r ) i.i.d. with p^(u r ), w r G {0, 1, • • ■ , [2 nR ' \ - 1} := S r 

. u[, 3) O r o) i.i.d. with p (3) (ii b ), w r o G 5 r0 

and define bin B { := {w r \w r G [(i - 1) • L2 n (^- Ra °)j + l,i • L2 n (^" Ra °)j]} for t G 5 a0 . 

Encoding: During phase 1 (resp. phase 2), the encoder of node a (resp. b) sends the codeword x. a (w a ) 
(2) 

(resp. x b (wb)). At the end of phase 1, relay r estimates (or decodes) w a . At the end of phase 2, the 

(2) (2) (2) 

relay compresses the received y r to a message w r o if there exists a w r o such that (y r , yy (io r o)) G 
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Fig. 13. Maximum sum-rate of the 8 bi-diretional protocols 
and 2 outer bounds at different SNR. Here h ar — 0.6, h br = 
20 and h ab = 0.5 



Fig. 14. Maximum sum-rate of the 8 bi-diretional protocols 
and 2 outer bounds at different SNR. Here h x = 20, h b , = 
0.6 and h 3b = 0.5 




Fig. 15. Maximum sum-rate of the 8 bi-diretional protocols at different relay position. Here ft a b = 0.2 and no rate constraints. 



(Y r Y r ). Such an w r o exists with high probability if 

ii r o = A 2in /(y r (2) ;f r (2) ) + e (147) 
and n is sufficiently large. Also we choose 

i? r = A 3 ,„/([/ r (3) ;Y b (3) )-4e (148) 
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Fig. 16. Maximum sum-rate of the 8 bi-diretional protocols 
at different relay position. Here h a b = 0.2 and a = 1 (R 3 — 
R b ). 



Fig. 17. Maximum sum-rate of the 8 bi-diretional protocols 
at different relay position. Here /i a b = 0.2 and a = 2 (_R a = 
2R b ). 



To choose w r , the relay first selects w 3 o = s a (w a ) and the bin B Wa0 , then it searches the minimum 
w r G B Wm such that (ii, (^r)i u h (^ro)) £ A®(l7 r l7b). This ensures uniqueness of u> r if such a u> r 
exists, i.e., u> r is a function of (w a o,w r o). Such a w r exists with high probability if 

| 5saK) | =2 -(A 3 ,„/(C/^ ; ^)+ e ) (149) 

Since |J3j| = 2 n ( fir_i?a °), Vj G <S a o> the condition is equivalent to 

R a0 < A 3 ,„/(C/ r {3) ; y b {3) ) - A 3 ,„/(C/ r (3) ; C/^ 3) ) - 5e (150) 

The relay then sends x^ 3 " 1 randomly generated i.i.d. according to p( 3 \x r \u r , ut>) with u[ 3 ^(u> r ) and 
u b C^ro) during phase 2. 

Decoding: Node a estimates ?D r o after phase 3 using jointly typical decoding. First, since node a 
knows w a , it can reduce the cardinality of w r to |^2 n (- R ' _ - Ra °) J . Furthermore, it forms two sets of w r o 
based on typical sequences, {«; r o|(ya ,yp (w r o)) G ^ (2) (^a^r)} and {wro|(up(u; r ), u b 3) (u} r0 ),yi 3) ) G 
^4^ 3 ^ (C/ r C/b^ a ) 5 f^r G S TOa0 }. After decoding zZ; r o (which is a success if there is a single common ele- 
ment in both of the previous sets), node a decodes using jointly typical decoding of the sequence 

(2) (2) (2) 

(x. , y r ' , yJj ). Node b decodes w r after phase 3 and from the bin index of w r it estimates s a . Then node 
b decodes the index as w a if there exists a unique w a G 5 a ,s a such that (x a ^ (w a ) , y ^ ) G ^4Wpf a Yb)- 
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Error analysis: 



P[E ijb ] < P[Ei]} U E® UE®} (151) 

< P[E$] + P[E^ \E$] + P[E® \e£) n E®] (152) 

P[£ b , a ] < P[45 U £g] (153) 

^P^g'l + P^gl^J] (154) 

We define error events in each phase as follows: 

1) eM=eM' 1 ue$> 2 . 

Eif 1 : (x«K),y«)^AW(X a y r ). 

E^)' 2 : there exists w a ^ w a such that (xjj 15 («;,), yf 1 *) € AW(jr a y r ). 

2) = i^' 1 U E^> 2 U £ r (3) ' 3 U 



Ef)' 1 : there does not exist a iu r o such that (y[ 2) ,y( 2) Ko)) G A( 2 )(r P Y r ). 



p(3),3 

r,b 
P (3),4 



there does not exist a u; r € B Wm such that (u[ 3 ^ (w r ), u b (w r o)) G A( 3 )([y r [/ b ). 
(ui 3) (w r ),y®)?A®(U r Y h ). 

there exists w r / w r such that (ur (w r ),y b ) G A^{U r Y b ). 



on P (3) _ p (3),l p (3),2 
3 ) ^a,b - ^a,b U ^a,b 



E^l' 2 : there exists w a ^ w a such that (x^(it; a ),y^) G ^W(X a y b ) and w a € S- 



4) Eg = ES' 1 U £$' 2 U 4 3) ' 3 U E^ A U Eg' 5 U £ r (3) < 6 . 

Eg' 1 : there does not exist a u; r0 such that (y^.y^Ko)) G A^(Y r Y r ). 



F (3),2 

p(3),3 
^r,a 

F (3),4 
^r,a 

F (3),5 



there does not exist a u; r € S Wa0 such that (iir (ii? r ), uj^ (itvo)) G ^4^(f7 r C/b)- 

(yi 2) ,y r (2) Ko))0A( 2 )(y a y r )- 
(u( 3) K),u b 3) K o ),yi 3) )0yi( 3 )(^^ b y a ). 

(3) (3) (3) 

there exists (w n w r o) where w r / iu r andw r o ^ uvo such that (uj (u> r ),u b '(w r o),y a ) 



A®(U r U b Y a ) and (yi 2 \y r (2) (w r o)) G A( 2 )(Y a Y r ). Recall, w r is uniquely specified by 
(w a o, Wro). Hence for a given w a o, there are at most 2 nRr0 such (w r ,w r o) pairs. 
E r (3) ' 6 : there exists w r0 + w r0 such that (u r (3) K), u b 3) (u} r0 ), y£ 3) ) G A( 3 )([/ r £/ b Y a ), 

(yF,y r (2) Ko)) GA( 2 )(y a y). 



5) E$=E^UEir. 



E ( b V ■ (x b 2 \w b ),y?\y?\w r0 )) ? A&(X b Y a Y r ). 
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P®' 2 : there exists w b / w b such that (x {2) (w b ),y {2) (w r0 )) E ^ 2 ) (X b Y a Y~ r ). 
Then, 

P[P$] <P[E^/] + P[P^' 2 ] (155) 

<P[5( 1) (x a K),y r )]+P[U (iaM Z)( 1) (x a (ti a ),y r )] (156) 

< e + 2 ™(^-Ai,„/(x< 1) ;y r < 1) )+3£) (157) 

P[E®\e£}\ ^PiE^lBQ] + P[E<g' 2 |£#] + P[E®' 3 \Ej8] + P[£ r ( J' 4 |£$] (158) 

<2e + P[5( 3 )(u r K),y b )] + P[U^ rJ D( 3 )(u r (u; r ),y b )] (159) 

< 3e + 2 n(R r -A3,„/(f/< 3 );y b < 3 ))+3 e ) (160) 

In (fl59l . PfP^' 1 !^] and P[£;JJ ),2 |.E$] are less than e due to (fl47T) and (fBOl . respectively. 

P[E$>\M$nM$] ^[^lE^n^+P^l^n^] (161) 

<P[P« (x a («;,), y b )] + P[U^ a ^ a P (1) (x a (w a ), y b ), s a (w a ) = s a (w a )] (162) 

< e + 2 ™(^-Ai,„/(X< 1) ;Y b (1) )- J R a0 +3e) (163) 

P[P r ( J] ^PfPg' 1 ] + PfPg' 2 ] + PfPg' 3 ] + P[P r (3 a )>4 ] + P[E$' 5 ] + P^' 6 ] (164) 

<2e + P[ J D^(y a ,y r Ko))]+P[ J D (3) (u r K),u b Ko),y a ))] + 
P[U Z?^(u r (u; r ),u b (u; r o),y a ),L> (2) (y a ,y r (u; r o))] + 

P[U «,=., Z)( 3 )(u r (u; r ),u b (u; r o),y a ), J D^(y a ,y r (u; r o))] (165) 

<3 e + 2 ™(fto-A3,„/(C/ r (3) ,!7 b <3) ;ya (3) )-A 2 ,„/(y a (2) ;K <2) )+7e) + 

2 n(iJ rf) -A 3 ,„/(^ b <3) ;C/ r <3) ,y a (3) )-A 2 ,„7(Y a (2) ;y r (2) )+7e) (166) 

In (f!65l PfP^' 1 ] and PP^* 2 ] are less than e due to (fl47T) and (fBUb . respectively. In (fl66l P[L>( 2 ) (y a , y r (u>ro))] 
is less than e by the Markov lemma. 

P[Pg 1^3] < P l E bT\ B ^ + ^^S'^rJ] (167) 

<P[D( 2 )(x b (u; b ),y a ,y r Ko))]+P[U^ b D( 2 )(x b (u; b ),y a ,y r K ))] (168) 

< e + 2 n(flb-A 2i „/(X b (2) ; y a (2 ',Y r (2) )+4 £ ) (169) 

Since e > is arbitrary, the proper choice of P a o, the conditions of Theorem [Q] (1 148b . (11501 ) and the 
AEP guarantee that the right hand sides of (1 1 57b . (II 601 ). (11631 ). (11661 ) and d 1 69b corresponding to the 
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first term of (|84ll . (11481 ). (the second term of dHJ and (11501 1 s ). (l86l ) and dHSb vanish as n — > oo. By the 
Caratheodory theorem in [9], it is sufficient to restrict \Q\ < 6. ■ 
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